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Measurements of production cross sections and polarizations are essential inputs for a theoret- 
ical understanding of heavy vector meson production. In this article the CDF measurement of 
the Y(\S) polarization in the kinematic range \y\ < 0.6 and 2 < pr [GeV/c] < 40 using a data 
sample of 2.9 fb _1 is described. Compared to the CDF Run I measurement, with which it agrees, 
it extends the covered pj range, allowing for a better comparison to predictions in the perturba- 
tive regime. The observed trend towards longitudinal polarization at high transverse momentum 
disagrees with predictions based on non-relativistic QCD. 
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1. Introduction 

The understanding of heavy vector meson production at hadron colliders and ep colliders has 
evolved in the past several years. While the Color Singlet Model describing the transition for 
the heavy quark anti-quark pair produced in the hard scattering process to the bound QQ state by 
the emission of hard gluons, failed to reproduce the J/y cross section measured at the Tevatron 
[g], the addition of soft gluon emission in the Color Octet Model [Q] could solve this issue. The 
theoretical framework usually applied here is non-relativistic QCD (NRQCD) It factorizes the 
quarkonium creation into a hard process to produce the QQ pair, and a long-distance matrix element 
for the meson formation. As the non-perturbative matrix elements are believed to be universal, they 
can be measured in experiments and then applied in theoretical predictions. 

One of the predictions of NRQCD is that the polarization of heavy vector mesons should be 
transverse at large transverse momenta. Recent CDF measurements of the J/\j/ polarization [^] 
disagree with this expectation. A possible explanation for this discrepancy is that the charm quark 
is too light for perturbative calculations. Fortunately, the Tevatron provides the possibility to test 
the predicted trend towards transverse polarization for the much heavier Y states 

2. Spin Alignment Measurement 

The polarization parameter a is defined as a = gp^f^ where Oj/l are the cross sections 
for the production of transverse and longitudinal polarized vector mesons, respectively. Thus the 
minimal and maximal values of a = — 1 and a = 1 correspond to full longitudinal and transverse 
polarization, respectively. A value of a = indicates unpolarized production. 

The measurement of the spin alignment of the vector meson requires the definition of a refer- 
ence axis. In this measurement the reference axis is defined by the Y momentum in the lab frame. 
This definition, called s-channel helicity frame, is used by all polarization measurements at hadron 
colliders so far. 

The polarization parameter a can be determined experimentally by measuring the distribution 
of the decay angle 6* of the positively charged muon from the Y — > pL + pi~ decay in the Y rest frame 
with respect to the reference axis. The angular distribution is given by d ^ e ^ x 1 + ot cos 2 6*. Be- 
cause cos 8* enters only quadratically, the distribution is symmetric about 0. This consequence of 
parity conservation allows to perform the measurement as function of | cos 6* | and thus to increase 
the statistics per bin in the angular distribution. 

3. Data Sample 

The analyzed data sample corresponds to an integrated luminosity of 2.9 fb _1 of pp collisions 
at a center of mass energy of y/s = 1 .96 TeV recorded by the CDF II detector at the Tevatron 
collider. The Y events are triggered by a pair of oppositely charged muons with a transverse mo- 
mentum of Pt(h) > 3 GeV/c. In the offline selection one muon is required in addition to have 
at least a transverse momentum of 4 GeV/c. Only muons in the central region of the detector, 
with a pseudorapidity |tj(/i)| <0.6, are used. Further selection criteria are applied to select well- 
reconstructed Y mesons. The kinematic range of Y mesons studied in this analysis covers the 
central rapitidy region, \y\ < 0.6, and the transverse momentum range 2 < pj [GeV/c] < 40. 
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Figure 1: Invariant jJ. + jJ.~ mass distribution with shaded Y(IS) signal and side band regions. 

Figure [I] shows the invariant mass spectrum of the selected dimuon pairs together with a fit 
where the signal shapes are taken from simulation and the background is described by an expo- 
nential. A range of ±2. 5a around the F(l.S') mass peak center is defined as signal region. The 
good mass resolution allows to define lower and upper mass sidebands close to the signal for the 
estimation of background events in the signal region based on the exponential background shape. 
The selected data sample contains about 80.000 Y(IS) mesons. 

4. Analysis Method and Systematic Uncertainties 

The same method to determine the polarization as used in Run I [[7|] and for the Run II y 
polarization measurement [^] is applied here. The data sample is split into bins of pj and then the 
cos0* distribution in each pj bin is examined. A % 2 fit is performed to the | cos 0* | distributions 
with 10 equidistant bins. The background contributions for all bins are free parameters in the 
fit and constrained by simultaneously fitting the mass sideband regions. The signal distribution is 
described by the weighted sum of two templates for transverse (T) and longitudinal (L) polarization. 
The relative weights of both contributions determine the polarization parameter a. 

Because the |cos0*| signal templates are affected by acceptance and efficiency effects, they 
are determined from simulation. Fully-transverse and fully-longitudinal polarized Y(IS) — > Ai + jU~ 
samples were generated with the EvtGen package and simulated and reconstructed with the stan- 
dard CDF software. Trigger efficiency effects were taken into account by applying efficiencies 
measured on data to the simulation. Since the acceptance causes the | cos 8* | template distributions 
to be dependent on the transverse momentum of the Y, the simulated pj distribution is reweighted 
in an iterative procedure to match the one observed in data. Figure || shows the fit for one pj bin 
and the obtained background distributions for low, medium, and high p T . The background shows 
a strong dependence on the transverse momentum. 

Several systematic effects that could alter the | cos 0* | distributions were investigated. The 
uncertainty on the trigger efficiency measured on data leads to a systematic uncertainty of 0.007 
on the polarization parameter. Other systematic uncertainties were found to be negligible. This 
includes variations of the mass shapes and the signal window definition, a change of the Pt{h) cut, 
and using an alternative pr reweighting. Although the chosen bin width is much larger than the 
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Figure 2: Fit (black) to |cos0*| distribution (green) for 2 < p T [GeV/c] < 3 with T (red) and L (blue) 
templates (left) and fitted angular distributions of background events for selected pj bins (right three plots). 
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Figure 3: Y(IS) polarization measurement compared with the CDF Run I result (left) and a NRQCD pre- 
diction M (right). 



resolution in cos (9*, a check for a systematic effect was performed by using twice the number of 
bins. No effect beyond expected statistical fluctuations was found. It was also verified that there is 
no sizable migration between pj bins. Overall the systematic uncertainties are much smaller than 
the statistical ones. 



5. Result 

Figure [| shows the measured Y(IS) polarization as a function of pj. It is compatible with 
zero polarization, with a trend towards longitudinal polarization at high pj. As can be seen in the 
left plot, the result is consistent with the CDF Run I measurement [f7|] and considerably extends its 
momentum range to higher pj. The right plot in Fig. |3| illustrates that the measurement contradicts 
the NRQCD prediction of transverse polarization at high pj. 

As shown in Fig. |4] the result disagrees with the DO measurement Since the DO measure- 
ment indicates a longitudinal polarization at low pj, the result in this region was cross-checked 
by comparing the measured |cosf3*| distribution with the simulation of unpolarized events. Good 
agreement is observed as shown exemplarily for one pj bin in the right plot in Fig. ||, confirming 
the consistency with zero polarization. 
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Figure 4: Y(1S) polarization measurement compared with the DO result (left) and comparison of data 
with simulated unpolarized events at low pj (right). 

A possible reason for the discrepancy between the CDF and DO results is that the DO mea- 
surement was performed in a larger kinematic range of \y\ < 1.7. Another difference between both 
measurements is caused by the different track reconstruction precisions. Due to a worse mass res- 
olution the signal shapes of Y(IS), Y(2S), and Y(3S) overlap in the DO measurement, not allowing 
to define mass sidebands as close to the signal as in the CDF measurement. 

6. Summary and Outlook 

The Y(IS) polarization is measured by CDF up to pj = 40 GeV/c. It is consistent with zero 
and tends towards longitudinal polarization at high pj. The result is compatible with the CDF run 
I result and inconsistent with the DO result. It also disagrees with NRQCD predictions. 

Updated measurements with more data, having data with about three time the integrated lu- 
minosity already recorded, can be expected, including polarization measurements of the Y (2S) and 
F(35). This, together with a separation of direct and feed-down contributions, may help to obtain 
a better theoretical understanding of heavy meson production. 
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